Near-optimal dynamical decoupling of a qubit 



Jacob R. West, 1 Bryan H. Fong, 1 and Daniel A. Lidar 2 

tfifRL Laboratories, LLC, 3011 Malibu Canyon Rd., Malibu, California 90265, USA 
2 Departments of Chemistry, Electrical Engineering, 
and Physics, Center for Quantum Lnformation & Technology, 
University of Southern California, Los Angeles, California 90089, USA 
(Dated: July 7, 2009) 

We present a near-optimal quantum dynamical decoupling scheme that eliminates general deco- 
herence of a qubit to order n using 0(n 2 ) pulses, an exponential decrease in pulses over all previous 
decoupling methods. Numerical simulations of a qubit coupled to a spin bath demonstrate the 
superior performance of the new pulse sequences. 
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Quantum information processing requires the faithful 
manipulation and preservation of quantum states. In the 
course of a quantum computation, uncontrolled coupling 
between a quantum system and its environment (or bath) 
may cause the system state to decohere and deviate from 
its desired evolution, potentially resulting in a compu- 
tational error. Here we present a dynamical decoupling 
(DD)[1 scheme designed to mitigate this effect. DD com- 
bats this decoherence by suppressing the system-bath in- 
teraction through stroboscopic pulsing of the system, an 
idea which can be traced to the spin-echo effect, with 
a long tradition in NMR [2]. Our new DD scheme is 
near-optimal, and provides an exponential improvement 
over all previously known DD protocols. It suppresses 
arbitrary coupling between a single qubit and its envi- 
ronment to n th order in a perturbative expansion of the 
total qubit-bath propagator, using 0(n 2 ) pulses. 

The most general interaction between a qubit and a 
bath can be modeled via a Hamiltonian of the form H = 
I^Bj+X^Bx+Y^By + Z^Bz, where X, Y and Z are 
the Pauli matrices, / denotes the identity matrix, and B a 
are arbitrary Hermitian bath operators. The term Bj is 
the internal bath Hamiltonian. Our scheme builds upon 
two recent insights. The first is due to Uhrig [3], who - 
for singlc-qubit decoherence consisting of pure dephasing 
errors (Hz = I "® Bq + Z (g) Bz) - found a scheme ("Uhrig 
DD" - UDD) which prescribes a sequence of X pulses 
(7r-rotations around the x axis) at times 
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where T is the total evolution time and j = 1, 2, . . . , n, if 
n is even, and j = 1, 2, . . . , n+1, if n is odd. These times 
characterize a filter function [4, 5 that removes the qubit- 
bath coupling to n th order in a perturbative expansion 
of the total system-bath propagator. UDD is provably 
optimal in that it achieves the minimum number of pulse 
intervals, n+1, required to accomplish this removal [HI [7]. 
For single-qubit decoherence consisting of pure spin-flip 
errors (Hx = I®Bq+X®Bx), UDD is still optimal, sim- 
ply by substituting Z pulses (7r-pulses around the z axis) 
for X pulses. However, for systems subject to general er- 
rors as prescribed by H , DD schemes incorporating both 



X and Z pulses are required. One such scheme, which 
provides our second source of insight, is concatenated 
DD (CDD) [5]: the CDD sequence is capable of elimi- 
nating arbitrary qubit-bath coupling to order n at a cost 
of 0(4 n ) pulses [9]. CDD works by recursively nesting 
a pulse sequence found in Ref. [T], capable of canceling 
arbitrary decoherence to first order. Uhrig recently intro- 
duced a hybrid scheme (CUDD) which reduces the pulse 
count to 0(n2 n ) for exact order n cancellation [TU]. By 
appropriately concatenating the UDD sequences for Hz 
and Hx we show here how arbitrary decoherence due to 
H can be exactly canceled to order n using only (n+1) 2 
pulse intervals. A numerical search we conducted found 
that this is very nearly optimal for small n, differing from 
the optimal solutions by no more than two pulses. 

Near-optimal pulse sequence construction. — The goal 
of our construction is to integrate an X-type UDD„ se- 
quence, which suppresses pure dephasing error to order n, 
with a Z-type UDD„ sequence, which suppresses longitu- 
dinal relaxation to order n, so that the resulting sequence 
removes arbitrary error to order n. Since the total time 
T in Eq. (|T|) is arbitrary, what matters for error cancel- 
lation is not the precise pulse times tj, but rather the 
relative sizes of the pulse intervals Tj = (tj — tj-i), for 
j = 1,2, ...,n + 1. Thus, the most relevant quantities 
are normalized pulse intervals, 
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again with j = 1, 2, . . . , n+1. Here we chose to normalize 
with respect to the shortest pulse interval (t\ — to) = t\, 
so that Tj = Sjt\, which has the important consequence 
that the (normalized) total time grows with n, as pulse 
intervals are added to address higher order error. This 
fixing of the minimum pulse interval ti corresponds to 
imposing a finite bandwidth constraint. Since any phys- 
ical implementation will not be able to shrink the pulse 
intervals arbitrarily, but will be limited by the fastest 
pulsing technology available, the change in perspective 
from fixed total time to fixed minimum interval is appro- 
priate. The total normalized time of a UDD„ sequence 
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is then given by, 
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so that the total physical time is T = SVA- 

Let U{t) denote the joint unitary free evolution of a 
qubit and its bath for a time r, subject to the Hamilto- 
nian H . A Z-type UDD„ sequence then takes the form, 

Z n (r) = Z n U(s n+1 T) Z U( Sn r) - Z U(s 2 r) Z U( Si t). 

(4) 

Note that a final Z pulse is required for n odd. 
Define the X-type UDD„ sequence X n (r) similarly. 
Yang and Liu showed in [B] that the entire UDD 
sequence X n (ti) can be expressed as the propa- 
gator exp (-iS n h(I ®B' I + X®B' X ) + OdXSnh)^ 1 )) 
(in units of h = 1), where A ~ ||-ff||, provided XS n ti 
is sufficiently small to ensure convergence of the time- 
perturbative expansion. The important point here is that 
the resulting effective Hamiltonian H' = I eg) B\ + X ® B' x 
is a purely spin-flip, time-independent, coupling. The 
correction term 0((XS n ti) n+1 ) potentially contains all 
couplings with a complicated time-dependence, but is 
suppressed to order n. Moreover, subject to the con- 
vergence condition, we are free to choose the minimum 
interval arbitrarily without impacting the validity of the 
proof in [5] . This is the key to correctly integrating the 
X-type and Z-type UDD„ sequences. The desired er- 
ror cancellation properties only require that the normal- 
ized pulse intervals sj have the specified structure. The 
precise physical timing is inconsequential. Therefore, to 
integrate X n and Z n sequences properly, without break- 
ing the delicate pulse timing structure required for error 
cancellation, we must scale the pulse intervals of the in- 
ner sequences uniformly with respect to the outer pulse 
sequence structure. Hence, if each U(sjr) in Eq. Q is 
replaced by the time-scaled DD sequence X n (sjT), then 
the outer Z n sequence suppresses the purely spin-flip cou- 
pling H 1 remaining after each X n {sjT) sequence, produc- 
ing general decoherence suppression to order n with only 
(n + l) 2 pulse intervals. So the combined, near-optimal 
"quadratic DD" (QDD) sequence, takes the form, 

QDD„(r) = Z n X n (s n+1 T) Z X n {s n r) ■ ■ ■ Z X n ( Sl r), 

(5) 

abbreviated by the notation QDD„(r) = Z„(X„(t)). 
Notice how the relative scales of the pulse intervals 
are preserved. In each of the inner X n (sjT) se- 
quences, the ratio between successive intervals remains 
(sk+iSjT)/(skSjT) = Sk+i/sk, and for the outer Z n 
sequence the ratio is (S n Sj+iT)/(S n SjT) — Sj+i/sj, 
thereby ensuring the error cancellation properties of 
each sequence are left intact. Of course, an equivalent 
QDD„ sequence may be constructed as X n {Z n {fj). 
Moreover, though the inner DD sequences must have 



an equal number of intervals, they need not be the 
same length as the outer sequence, but can instead 
be adjusted to more efficiently address the dominant 
sources of error in any particular implementation. 
This way, QDD m n = Z m (X n (T)) is the more general 
construction, where the inner sequences suppress one 
type of error to order n, while the outer sequence 
suppresses the remaining error to order m. As the sim- 
plest explicit example, QDD X = ZX 1 (s 2 t)ZXi(sit) = 
Z (X £7(s|t) X U{ Si s 2 t)) Z {X U(s 2Si t) X U{s\t)) 
Y U(t) X U(t)Y U(t) X U(t), which we recognize as 
the so-called "universal decoupler" sequence found in pQ 
and used as the basis of the CDD sequence in [8] . 

It is worth noting that the construction of QDD„ af- 
fords a nice visualization. For a 6 {I, X,Y, Z}, define 
a(sj) = aU(sjT) a, then consider, 
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with the final QDD„ pulse sequence formed by reading off 
rows of the inner square from top to bottom, or columns 
from left to right. In other words, our construction may 
be succinctly described as an outer product between X- 
type and Z-type UDD sequences. 

Also, note that the QDD„ sequence requires a total 
physical time of S^T, and just as in the original proof 
[6] of order n error suppression for the Z n sequence, the 
condition for convergence of the perturbative expansion 
still remains, namely that XS^t is sufficiently small. An 
immediate consequence of this constraint is that there 
exists a maximal order of error n that can be suppressed 
for a given r, beyond which the error cancellation prop- 
erties of the sequence begin to break down. Conversely, 
if one hopes to suppress error to some fixed order n, then 
this implies a maximum pulse rate r = 1/r which must 
be attained. Specifically, a sufficient condition for con- 
vergence is XS^t < 1, hence r > Acsc 4 (7r/(2n + 2)). 

Finally, we translate these results back into a precise 
physical timing for the individual pulses. If the total time 
for the sequence is T, then the outer Z n sequence requires 
that Z pulses occur at the original Uhrig times prescribed 
in Eq. , while the inner X n sequences require X pulses 
executed at the times, tj_k = tj sin 2 ^ 2 n+? ) is where 
j,k = 1, 2, ... Ti if n is even and j, k = 1, 2, . . . 7i+ 1 if n is 
odd. When n is odd, X and Z pulses coincide at times 
tj, in which case Y — ZX pulses are used. 

Numerical results. — We now present numerical results 
that illustrate the efficiency of our new DD pulse se- 
quences in preserving arbitrary initial quantum states. In 
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the full-state quantum memory simulations that follow, 
the system and bath are initialized together as a (gen- 
erally nonseparable) random pure state The results 
presented involve a single system qubit coupled to four 
bath qubits. While this is clearly an unrealistically small 
bath, larger simulations we performed indicate that the 
relevant error suppression properties of our DD sequences 
are qualitatively unaffected by bath size. Similarly ideal- 
ized are the DD pulses themselves, which we take to be 
infinitely strong zero- width pulses, consistent with the 
analysis in [Fj . The evolution of the coupled system and 
bath in our simulations is governed by the generic Hamil- 
tonian H, with additional parameters {J, /?} included to 
control the coupling and bath strengths, respectively: 

H = 13 {I <g> Bx) + J (X ® B x + Y <g> B Y + Z ® B z ) , 

(6) 

by B a = 
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ct* <» crj.) , with a, k, I € {I,X,Y,Z}, i,j 
indexing the bath qubits, and randomly chosen coeffi- 
cients rjy 6 [0,1]. Note that these bath operators in- 
clude all 1- and 2-body terms, so that the system-bath 
Hamiltonian includes 2- and 3-body terms. The numeri- 
cal simulations are run with more than 100 digits of pre- 
cision and results are averaged over 10 random realiza- 
tions, with each instance randomly generating new bath 
operators and a new initial state. The vertical axes in 
each of these plots quantifies the DD sequence perfor- 
mance as log 10 (D), where D is the standard trace-norm 
distance | ||p(T) — p(0)||j between the evolved system 
state p(T) = tr B (\4>) (<f>\), with \cf>) = QDD„(r) \ip), and 
the initial system state p(0) = trsdV') (V'l)- Here tr^ 
is the partial trace over the bath. This distance mea- 
sure bounds the usual Uhlman fidelity F from above and 
below, 1 — D < F < y/l — D 2 , and is itself bound by 
the norm of a perturbative expansion of the propagator 
|llj . Indeed, in our simulations, D always goes like the 
norm of leading order term, as one might expect. Error 
bars are included at every point indicating the maximum 
deviation from the average distance measured. 

Figure [l] shows the performance of various DD schemes 
plotted against the number of pulses each requires. 
PDD„ is the n-times repeated universal decoupling se- 
quence found in [I], requiring 4n pulses (this sequence 
is capable only of first order decoupling); CDD„ is the 
n-times concatenated universal decoupling sequence, re- 
quiring 4™ pulses [HJ; CUDD„ is concatenated Uhrig DD, 
combining an n th order A-type Uhrig sequence with n 
concatenation levels of Z pulses, requiring n2™ pulses 
|10j . The coupling parameters J and (3 are fixed so that 
Jt = (3t = 10 -6 , given a shortest pulse interval t. Re- 
call that since the shortest pulse interval is held constant, 
the total evolution time grows with increasing n. Specif- 
ically, for QDD n the total evolution time is S 2 t. The 
evolution times of the other DD sequences also increase 
with n, but at different rates depending on how the num- 
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FIG. 1: Performance comparison of various DD schemes vs. 
number of pulses required. QDD n dramatically outperforms 
all known alternatives in both state preservation and number 
of pulses. 



ber of required pulse intervals scales with n. Notice the 
two most visible consequences of increasing total time 
in Figure [lj (1) at n = the evolution time is only 
SqT = r, so the dominant contribution to the overall fi- 
delity is Jt = 10~ 6 , which explains why the undecouplcd 
evolution still achieves excellent agreement with the ini- 
tial state in this plot; and (2) as the total time increases 
so does decoherence, explaining why the performance of 
PDD„ actually degrades with large n. 

Figure [2] shows QDD n performance as a function of J, 
the system-bath coupling strength. As JS 2 t approaches 



QDD„ : Distance vs. coupling strength 
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FIG. 2: QDD n performance dependence on coupling strength 
J, with fixed (3t = 10" 6 . 

1, the qubit decoheres so rapidly that DD has essentially 
no effect, and the distance between the initial and fi- 
nal states approaches its maximum of 1. Again, the 
n = line represents undecoupled free evolution, for 
which JS 2 t = 1 when Jt = 1, corresponding to the 
point in this plot where the log 10 (D) becomes zero and 
subsequently stays there. State preservation improves 
as n increases and another order of error is suppressed, 
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with the effect magnified as Jt decreases, though increas- 
ing total time counteracts the overall performance gain. 
This is evidenced in how, at each fixed Jt, the magni- 
tude of the improvement decreases as n increases, or in 
other words, in how the gap between lines narrows as n 
increases. The dotted line in both Figures 2j3 indicates 
when Jt = (3t. Across this transition point, the slope 
increases as the leading order error term changes from (3 
to J dominated in Figure [2j and vice versa in Figure [3] 
Indeed, when J < (3, the leading order term contributes 

which gives a slope of 1 as 



to distance as 



(n+l)\ 



a function of J for every n in Figure [2] On the other 
hand, when J > j3, the dominant contribution to perfor- 
mance becomes ^ J ("J_\y_ — > explaining the observed slopes 
of n + 1 in this plot. 

Figure [3] shows the performance dependence of QDD„ 
vs. P, the parameter which quantifies the pure bath en- 
ergy scale. We again see the Jt = [3t transition in this 



QDD„ : Distance vs. bath energy 
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FIG. 3: QDD n performance dependence on bath strength, 
with fixed Jt = 10" 6 . 



Conclusions. — We have presented a DD pulse se- 
quence construction that guarantees simultaneous can- 
cellation of both transverse dephasing and longitudinal 
relaxation to order n using only (n + l) 2 pulse inter- 
vals, or 0(n 2 ) pulses, an exponential improvement over 
all previous DD methods. The constructed sequences are 
near-optimal, in that an exhaustive numerical search we 
performed produced optimal pulse sequences of lengths 
7, 14, and 23, for n = 2,3, and 4, respectively, i.e., 
(n + l) 2 — 2 pulses. In view of this, we believe our con- 
struction is very nearly optimal, in addition to having 
the significant conceptual benefit of a simple algorithmic 
description. Moreover, it is worth noting that the fil- 
ter function may alternatively be optimized for specific 
bath noise spectra (Locally Optimized Dynamical De- 
coupling, LODD)[4j or over a range of bath frequencies 
(Optimized Noise Filtration through Dynamical Decou- 
pling, OFDD)[5 , resulting in pulse timings different from 
Eq. (JT|) and decoupling performance superior to both 
CPMG and UDD for high-frequency baths 010. Our 
construction can be easily applied to LODD and OFDD, 
but for ease of explanation we restricted our discussion to 
the pulse timing in Eq. . Finally, of particular impor- 
tance will be generalizing this work to incorporate finite 
width pulses and, ultimately, computation [12]. We look 
forward to experimental tests of these sequences. 
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figure (the dotted line), as well as the improvement with 
increasing n. The obviously significant feature of this 
plot however is the deviation, and then eventual return- 
ing, of the evolved state relative to the initial state as 
/3r sweeps through the regimes: (1) (3t < Jt < S~ 2 , 
(2) Jt < (3t < S~ 2 , and (3) (3t > S~ 2 . This behavior 
can be understood mathematically by considering an in- 
teraction picture with respect to the pure bath dynamics, 
governed by Hb = (3(1 ® Bj), then expanding the other 
bath operators in the Bj eigenbasis. Using a Dyson series 
expansion with JS 2 t < 1, we find that when (3t < Jt, 

(JS 2 t)' 1+1 

the distance goes like (n+\)\ — > as ^ n ^ ne P rev i° us plot, 
independent of (3; when Jt < (3t < S~ 2 , the distance 
goes like {jlJ ® (n+\)\ > gi vm S a sl°P e of n as a func- 
tion of (3; and when [3t > S~ 2 , the distance goes like 



(C is a constant), which gives a 



I) < 

p 1 1 > (n+l)! 

slope of —1 until (3 is sufficiently large as to eliminate 
the contribution of the first term, leaving a leading order 
term that goes again like the small (3 regime. 



[1] L. Viola, E. Knill and S. Lloyd, Phys. Rev. Lett. 82, 2417 
(1999). 

[2] C. Slichter, Principles of Magnetic Resonance, Springer- 
Series in Solid-State Sciences (Springer, Berlin, 1996). 

[3] G.S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 

[4] M. J. Biercuk et al, Nature 458, 996 (2009). 

[5] H. Uys, M. J. Biercuk, and J. J. Bollinger, Phys. Rev. 
Lett. 103, 040501 (2009). 

[6] W. Yang and R.-B. Liu, Phys. Rev. Lett. 101, 180403 
(2008). 

[7] Optimality follows from the fact that there are n alge- 
braic constraints satisfied by n + 1 pulse intervals for 
0(n) decoupling. 

[8] K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 
180501 (2005). 

[9] K. Khodjasteh and D. A. Lidar, Phys. Rev. A 75, 062310 
(2007). 

[10] G.S. Uhrig, Phys. Rev. Lett. 102, 120502 (2009). 

[11] D. A. Lidar, P. Zanardi, and K. Khodjasteh, Phys. Rev. 

A 78, 012308 (2008). 
[12] K. Khodjasteh and L. Viola, Phys. Rev. Lett. 102, 



080501 (2009). 



